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In this paper we study the quantum properties of the lattice vacancy in topological band insulators 
(TBIs) and topological superconductors (TSCs) on honeycomb lattice with particle-hole symmetry. 
Each vacancy has one zero mode for the Haldane model and two zero modes for the Kane-Mele 
model. In addition, in TSCs on honeycomb lattice with particle-hole symmetry, we found the 
existence of the Majorana zero modes around the vacancies. These zero energy modes are protected 
by particle-hole symmetry of these topological sates. 



Topological band insulators (TBIs) represent a class 
of novel states of matter characterized by a special band 
structure : they are the materials that behave like insula- 
tors in their interior or bulk while permitting the metallic 
boundaries In two dimensions, there are two types 
of TBIs, the TBIs with time-reversal symmetry (TRS) 
and those without TRS. For the TBIs without TRSjl], 
the Thouless-Kohmoto-Nightingale-Nijs (TKNN) num- 
ber can be identified to be the topological invariant Q; 
While for the TBIs with TRS (quantum spin Hall states 
in two dimension[3, Q), people use the Z2 topological in- 
variant to label their properties Q- On the other hand, 
recently, people found that superconducting (SC) states 
with the same local order parameter may have different 
topological properties, which leading to the concept of 
"topological superconductivity (TSC)" d [T^]. 

The non-trivial topology of the bulk can also be ex- 
posed by introducing topological defects, tt-Aux, dislo- 
cations and vortices, to name a few. Moreover, the 
interplay of defect topology and topology of the origi- 
nal states result in even richer phases. Previous works 
have been done on the effect of tt-Aux and dislocations 
on the Haldane model which indeed trap zero energy 
bound stat es flTi IT^. Moreover, according to the index 
theorem jl3l4l5| . a vortex in topological superconductors 
is predicted to harness a Majorana fermion zero mode 
inside the vortex coreQ. In turn, such vortices were 
found to be non-Abelian anyons which are the holy grail 
of topological quantum computation [loj . 

In this paper, we will focus on the effects of different 
type of defects-vacancies, in two dimensional topological 
states with particle-hole symmetry on the honeycomb lat- 
tice. FIG.l.a is the illustration of a vacancy on the hon- 
eycomb lattice. It is known in graphene, zero modes of 
electrons were found around the vacancies [l^ Here 
we also found that the lattice defects of the TBIs on 
honeycomb lattice have nontrivial quantum properties : 
for the Haldane model each vacancy has one zero mode 
and for the Kane-Mele model each vacancy has two zero 
modes. These zero energy modes in TBIs are protected 



FIG. 1: (Color online) (a) The lattice defect of the honeycomb 
lattice - vacancy, (b) The particle density of the zero mode 
around the vacancy of the Haldane model with particle-hole 
symmetry (e = 0). 



by particle-hole symmetry and the finite energy gap of 
the electrons. For the graphene system without energy 
gap. the zero modes (localized states) of electrons around 
the vacancies arc fragile. By adding the perturbations 
that break the particle-hole symmetry (for example, the 
staggered potential), these zero modes (localized states) 
disappear. On the other hand, our results are much dif- 
ferent from those in the Ref . [l^l , in which the effect of lat- 
tice vacancies in graphene with intrinsic spin-orbit inter- 
action is studied in a continuum model. And these zero 
modes around a vacancy cannot be simply regarded as 
the remnant of the gapless edge states in a continuum ef- 
fective model as people have done in Ref. [l^ . In addition, 
in TSCs on honeycomb lattice with particle-hole sym- 
metry, we found that there also exists (non-topological) 
Majorana zero modes around the vacancies. 

Vacancies for the Haldane model: We start from the 
spinless Haldane model, of which the Hamiltonian is de- 
fined as 



;cj + h.c. 



-t' Y: e^'-clc, 
-^^4ci-£Yclci (1) 
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where t and t' are the nearest-neighbor hopping 
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and the next-nearest-neighbor hopping, respectively. 
gi0ij jg complex phase of the next-nearest-neighbor 
hopping, and we set the direction of the positive 
phase to be clockwise = f)- £ is the co- 

efficient of the on-site staggered energy. Using the 
Fourier transformations, we can get the spectrum of 

free fermions as Ek = i-y/ (^^ — e)^ + |^fc|^ where 

ICfel = + 2 cos (VSfc.y) + 4cos cos 

and = 2t' [-2cos(3fc^/2)sin(^/3fc,y/2) -I- sin(v^fcy)] . 
For this free fermionic system, there are two phases: the 
TBI state with TKNN number C = ±1 and the nor- 
mal band insulator (NI) state. The phase boundary is 
2e = 6V3i'. 

Firstly we study the particle-hole symmetry of the 
Haldane model. After doing the particle-hole transfor- 
mation, c|g^ O — CigA and c\^g O Ci^B together with 
a complex (or charge) conjugate transformation (20i . the 
hopping terms are invariant but the staggered potential 
term is changed. So for the case of e = 0, the Hamiltonian 
have the particle- hole symmetry as = P'^HuP where 
P is the particle-hole transformation operator for spin- 
less fermion on the honeycomb lattice. Now we have the 
spectrum symmetry : each energy level of the electrons 
with positive energy E must be paired with an energy 
level of the electrons with negative energy —E. While 
for the case of £ ^ 0, the particle-hole symmetry is bro- 
ken, iJn 7^ P^HnP. However, now, we still have the 
spectrum symmetry a.s E —E. This is because after 
doing the particle-hole symmetry, the staggered poten- 
tial changes sign. Due to the translation symmetry, we 
may do a sublatticc transformation as A < — > B. Then 
the total Hamiltonian is invariant. 

We study the effect of a vacancy on the electronic 
states. It is obvious that such a defect that breaks trans- 
lation symmetry is a local distortion of the system and 
has no topological properties. We can consider a vacancy 
as a " hole" in the system by removing a site. Due to the 
non-zero TKNN number, there may exist topologically 
protected edge states on the boundary of this " hole" . It 
is known that at low energy the dispersion of edge states 
has a form as E{k) vpk where vp is Fermi velocity of 
edge states. When the size of the hole shrinks, the energy 
levels of the edge states become discrete and eventually 
the edge states on the boundary around the "hole" turn 
into localized states around the vacancy. Because the 
Haldane model is based on a bipartite lattice, when we 
remove a site to create a vacancy, there will exist an un- 
paired electronic state. Due to the spectrum symmetry 
{E —E), the corresponding localized unpaired elec- 
tronic state must have exact zero energy. Such zero mode 
is protected by the particle-hole symmetry of the TBIs. 

We calculate the electronic states for the Haldane 
model with a vacancy for the particle-hole symmetry case 
(e — 0) numerically. On a 36 x 36 lattice, we found a zero 
mode near the vacancy in the Haldane model. In FIG.l.b, 
we plot the particle density of the zero modes. The parti- 
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FIG. 2: (Color online) (a) The energy levels withthe spectrum 
symmetry (E <t4> —E): the localized state has zero energy 
shown by the blue line and the continuum, spectrum is shown 
by the blue region (b) The scheme of DOS of electronic states 
with a vacancy. In non-interacting limit, the contribution of 
DOS from the vacancy is a 5 function at a; = 0. Here we set 
t' = 0.1. 



cle density is localized around the defect center within a 
length-scale ~ {Af)~^ where A/ is the fermion's energy 
gap. FIG. 2. a shows the energy levels via t' for the Hal- 
dane model with a vacancy. Now the total number of the 
electronic states for this case is odd. FIG.2.b shows the 
density-of-state (DOS) of the electronic states, of which 
a S function is around u! = due to the zero mode of the 
vacancy. 

Besides, we calculate the induced quantum number 
on a vacancy. From the numerical results, we have 
found a zero mode on the vacancy for the Haldane 
model. So a vacancy possesses two localized states 
which are denoted by |-|-) and |— ) that are the oc- 
cupied and unoccupied state of fermions. respectively. 
Around a vacancy, the fermionic operators are expanded 
as V'(r, t) = Efc^o c^.e-^^^**, (r) + ^fe^o dle^''^'^l{r) + 
co^o{r), where Ck and dj, are the operators of fc 7^ 
modes that are irrelevant to low energy physics. ^'o('') is 
the wave-function of the unpaired zero mode, cq is the 
annihilation operator of the zero mode. We define the 
induced fermion number operators as Np = J : tp''tp : 
(Pv = cjco -I- X] (CfeCfc — d\dk) — \ where : 'ip'^il: : means 

normal product of tp''tfj. So we obtain the eigenvalues of 
induced fermion number operator as Np\±) = ±i|±). 
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FIG. 3: The energy of the locaUzed states around the va- 
cancies via the staggered potential e. For £ 7^ 0, there is no 
spectrum symmetry (E <S> —E) and the localized state has a 
finite energy shown by the blue line (the continuum spectrum 
is shown by the blue region). For e > 3-\/3i', the ground state 
is NI phase, there is no localized state around the vacancy. 



FIG. 4: The energy of the localized state around a vacancy 
via the on-site disordered parameter a. For a 7^ 0, there is no 
spectrum symmetry (E —E) and the localized state has a 
finite energy shown by the blue line (the continuum spectrum 
is shown by the blue region). 



The occupation (or unoccupation) of this zero mode lead 
to = f (or Np = — |) fractional electronic charge. 
Although there is no topological background, an induced 
fractional charge is trapped around each vacancy in the 
Haldane model. 

In addition, we found the parity effect for the Haldane 
model with multi- vacancy: for odd number vacancies (for 
example, 3 vacancies, 5 vacancies...), there exist zero en- 
ergy modes; while for even number multi-vacancy (for 
example, 2 vacancies, 4 vacancies...), the localized states 
will have finite energy. It is also the particle-hole sym- 
metry that enforces a zero energy level together with the 
pairs of energy levels with finite energy for the case of odd 
number vacancies. For example, for two separated vacan- 
cies, there exist two localized states (quasi-zero-modes) 
and the localized states slightly split due to the quan- 
tum tunneling effect between them. It is obvious that 
the energy splitting decays exponentially as the distance 
between two vacancies increases. 

We have found the particle-hole symmetry protected 
zero mode around vacancies in the TBI phase of the 
Haldane model. Then we add the staggered on-site po- 
tential which breaks particle-hole symmetry. When the 
particle-hole symmetry is broken, the localized state will 
shift from zero to a finite value. We plot the bound 
state energy via t' for the case of e = O.lt in FIG. 3. 
For e < Ec = 3^/3t' , in TBI phase, there exist a local- 
ized state around the vacancy. For e > Sc = 3\/3i', the 
ground state is NI phase with trivial quantum properties, 
then the localized states disappear. 

Finally, by varying local disordered potential on a 
given site iq gradually, the wave-function evolves from 
the extended state to the localized state. The Hamilto- 



nian model in Eq.(l) turns into 

-atY^ (c\c,+h.c)-at' e'^^-clc, 

+ Kjcj^g^Cjog^. (2) 

where Vq — {1~ a)t/ a. a is a parameter to tune the local 
disordered potential. In the limit of a — > 0, the on-site 
potential turns into infinite and the hopping parameters 
turn to zero around the site iq. On the other hand, in 
the limit of q — > 1, the system has translation invariance, 
of which the Hamiltonian reduces into Eq.(l) of £ = 0. 
By tuning a, we don't find any quantum phase transition. 
One extended state will turn into the localized state when 
there exists arbitrary small disordered potential with a 7^ 
0. See FIG.4. 

On the other hand, we studied the properties of the 
topological defects (the 7r-flux on a plaquette of the hon- 
eycomb lattice) in the Haldane model and also found a 
zero mode around it in the TBI phase. The zero modes 
around the 7r-flux are protected by the topological in- 
variable. In TBI phase, the localized state has exact zero 
energy which is robust against arbitrary perturbation. 
There also exists an induced fractional charge Np = § 
around the tt-Aux. With the induced fractional charge, 
each TT-flux carries a statistical angle 7r/4. 

Vacancies for the Kane-Mele model: The Hamiltonian 
of the Kane-Mele model is similar to that of the Haldane 
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model as 

HKM = -t (ct ,c,,, + /i.e.) - i' ^ e*^" cta,c, 

+ ^ H ^tC^.'T - ^ XI 4'tC*,<t (3) 
ieA.a- i£B,a 

where is the Pauh matrix and a are the spin-indices 
representing spin-up (a =t) and spin-down (<t =|) for 
electrons. Using similar calculations, we found that there 
exist two zero modes around a vacancy for the Kane- 
Mele model with e = and there also exists the parity 
effect for multi- vacancy case. These zero modes are also 
protected by the particle-hole symmetry. When there 
exist particle-hole breaking terms, the zero modes turn 
into the localized modes with finite energy. 

Then we calculate the induced quantum number on 
a vacancy for the Kane-Mele model. Using similar ap- 
proach, we found that for the TBI phase, the two zero 
modes around a vacancy correspond to four degener- 
ate energy levels denoted by |t-i-)® W-)® \\r-)i 
W-)® W+)® |i-)|2l|. Here | a+) and | (t_) are the 
occupied and unoccupied state, respectively. For the half 
filling case, the localized states around the vacancy are 
14,+), |t+)® li- )j of which there exists a spin-i 
moment as 

It-)® i;+) = \ It-)® i;+), 

lt+)® It-) = -\ lt+)® It-). (4) 

In Ref.jiS], people proposed that the Kane-Mele model 
can be realized in a Silicene based material with hon- 
eycomb lattice and strong spin-orbital coupling interac- 
tion. Then in this system people can detect the local spin 
moment by observing uniform spin susceptibility which 
obeys Curie- Weiss law as Xs ^ Nis{kBT)~^ where N\s is 
the vacancy number. 

Recently, people have studied the quantum properties 
of the TT-flux in Ref . [2l| - |23j . They found there exist two 
zero modes and induced spin-i moment on the 7r-flux. 

Vacancies for the TSCs on honeycomb lattice with 
particle-hole symmetry: In this part we will consider the 
effect of a vacancy in a C = ±1 chiral TSCs on honey- 
comb lattice with particle-hole symmetry. The effective 
model of the TSC is 

^^TSC = -ff H + AinducG X CiCj + k.C. (5) 

where Hn is the Hamiltonian of the Haldane model with 
e = and Ajnduco is the induced SC order parameter 



due to the proximity effect to an extended s-wave SC 
order. Now the ground state is really a C = ±1 TSC, 
of which the topological properties are similar to those 
of Px + ipy wave pairing TSC. Above effective model has 
particle-hole symmetry as, E —E. 

We then studied the properties of a vacancy by solving 
the Bogolubov-de Gennes (BdG) equations. On each va- 
cancy, we found a single zero-mode which is described by 
a real fermion field 7^ = / dr[uoip'^ + vqi/j] (7^ = 7) ob- 
tained as a solution of the BdG equations. When two va- 
cancies fuse (taken to the nearby sites in the honeycomb 
lattice), the result contains more than one quasi-particle 
due to the Ising fusion rule [25|: axa = I+ip where a de- 
notes a vacancy, / the vacuum and tJj the fermion. This 
fusion rule is same to that of the non-Abelian anyons. 
However, although trapping a Majorana zero mode, the 
vacancy has trivial quantum statistics and cannot be re- 
garded as a non-Abelian anyon. That is we found a sym- 
metry protected Majorana zero modes in two dimensions 
without topological defects. So we call it non-topological 
Majorana mode. 

Furthermore, we studied the properties of the tt-Aux for 
this system and also found a zero mode around of it in the 
TSC phase. So in the TSC state, the localized state has 
exact zero energy which is also robust against arbitrary 
perturbation. With trapping the Majorana zero modes, 
each TT-flux obeys the non-Abelian statistics. 

Conclusion: In the end we draw a conclusion. In gen- 
eral, for a TBI or a TSC on a bipartite lattice with 
particle-hole symmetry, there must exist zero energy 
modes around a vacancy and the parity effect for multi- 
vacancy. In this paper we investigate the lattice vacancy 
by using TBIs and TSCs on honeycomb lattice as ex- 
amples to show this effect. For the TBIs and TSCs on 
a TT-flux lattice (another two dimensional bipartite lat- 
tice), the properties of vacancies are similar to those 
on honeycomb lattice. As the remnant of the gapless 
edge state around a "hole" with minimize size, the lo- 
calized states around the vacancies will have exact zero 
energy which are protected by the particle-hole symmetry 
of these topological states. The perturbations breaking 
particle-hole symmetry will shift the zero energy of the 
localized state to a finite value. But the induced quantum 
numbers on the vacancies doesn't change until the energy 
gap closes and the quantum transition phase occurs. 

Finally we give a table to show the difference be- 
tween the zero modes around a topological defect (the 
TT-flux) and the the zero modes around a vacancy (non- 
topological defect) in difference topological states. 
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The Haldane model 


The Kane-Mele model The TSC on honeycomb lattice 


TT-flux 

vacancy 


Abclian anyon with statistical anf 
c/2 charge 


;le 7r/4 spin moment Non-Abelian anyon 

spin moment Majorana fermion mode | 



TABLE I: In different topological states, the yr-flux (topological defect) and the vacancy (non-topological defect) have different 
properties. 



K. V. Klitzing, G. Dorda, and M. Pepper, Phys. Rev. 
Lett. 45, 494 (1980). 

R. Prange and S. Girvin, The Quantum Hall Effect 
(Springer, New York, 1987); H. Aoki, Rep. Progr. Phys. 
50 (1987) 655. 

M. Z. Hasan and C. L. Kane, Rev. Mod. Phys. 82, 3045 

(2010) . 

X. L. Qi and S. C. Zhang, Rev. Mod. Phys. 83, 1057 

(2011) . 

F. D. M. Haldane, Phys. Rev. Lett. 61, 2015 (1988). 

D. J. Thouless, M. Kohmoto, M. P. Nightingale and M. 
den Nijs, Phys. Rev. Lett. 49, 405 (1982). 
C. L. Kane and E. J. Mele, Phys. Rev. Lett. 95, 146802 
(2005); 95, 226801 (2005). 

B. A. Bernevig, T. L. Huge and S. C. Zhang, Science 
314, 1757 (2006). 

G. E. Volovik, Zh. Eksp. Teor. Fiz. 94, 123 (1988) [Sov. 
Phys. JETP 67, 1804 (1988)]. 

N. Read and D. Green, Phys. Rev. B 61, 10267 (2000). 

C. Weeks, G. Rosenberg, B. Seradjeh and M. Franz, Na- 
ture Physics 3, 796 (2007); G. Rosenberg, B. Seradjeh, 
C. Weeks, M. Franz, Phys. Rev. B 79, 205102 (2009). 
V. Juricic, A. Mesaros, R.-J. Slager, and J. Zaanen, Phys. 
Rev. Lett. 108, 106403 (2012). 

S. Tewari, S. Das Sarma, and Dung-Hai Lee, Phys. Rev. 



Lett. 99, 037001 (2007). 
[14] L. Santos, Y. Nishida, C. Chamon, and C. Mudry, Phys. 

Rev. B 83, 104522 (2011 ). 
[15] R. Rov. larXiv:1001.2571l 

[16] Vitor M. Pereira, F. Guinea, J. M. B. Lopes dos Santos, 
N. M. R. Peres, and A. H. Castro Neto, Phys. Rev. Lett. 
96, 036801 (2006). 

[17] M. M. Ugeda, I. Brihuega, F. Guinea, J. M. Gomez- 
Rodriguez, Phys. Rev. Lett. 104, 096804 (2010). 

[18] M. Inglot and V. K. Dugaev, J. Appl. Phys. 109, 123709 
(2011). 

[19] W. Y. Shan, J. Lu, H. Z. Lu, and S. Q. Shen, Phys. Rev. 

B 84, 035307 (2011). 
[20] C.N. Yang and S.C. Zhang, Mod. Phys. Lett. B 4, 759 

(1990) . 

[21] S. P. Kou, Phys. Rev. B 78, 233104 (2008). 
[22] X. L. Qi and S. C. Zhang, Phys. Rev. Lett. 101, 086802 
(2008). 

[23] Y. Ran, A. vishwanath and D. H. Lee, Phys. Rev. Lett. 

101, 086801 (2008). 
[24] C. C. Liu, W. Feng, and Y. Yao, Phys. Rev. Lett. 107, 

076802 (2011). 

[25] G. Moore, and N. Read, Nucl. Phys. B 360, 362-396 

(1991) . 



